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^^ , We discuss the construction of perfect fluid stellar objects having optical geometries with 

multiple necks corresponding to spatially closed unstable lightlike geodesies. We prove that 

Cm I there exist physically reasonable models with arbitrarily many necks. We also show how a 

first order phase transition can give rise to quite pronounced secondary double necks. The 

.J analysis is carried out using a modification of a recent dynamical systems formulation of the 

K^ , TOV equations due to Nilsson and Uggla. Our reformulation allows for a very general family 

^•f^ • of equations of state including, for example, phase transitions. 

o 
(^ 

O ' 1 Introduction 

o' 

^^ ' The theory of static spherically symmetric (SSS) relativistic stellar models has received an increasing interest 

qh, in recent years. A few years ago Chandrasekhar and Ferrari pi discovered the phenomenon of trapped 

gravitational waves in ultracompact objects. Later Abramowicz et al. [E| showed that the trapping could 
t^J^, be conveniently visualized in terms of the optical geometry of the models. It follows from their work that 

trapping occurs when the optical geometry has a sufficiently pronounced neck. The physical relevance of 
the neck stems from the fact that it is associated with an unstable closed lightlike circular orbit. This 

/^ ' phenomenon is present in the exterior Schwarzschild geometry at i? = iM as is well-known. Any SSS stellar 

object which is more compact than R = 3M will therefore have a neck. However, Rosquist showed in Q 
that the limit R < 3M is not essential for the formation of a neck and hence it is not essential for a trapping 
region either. In particular, a neck can exist in the stellar interior. Because the optical geometry is regular 
all the way down to the center of the star, a neck is always associated with a bulge which is located further 
inside the stellar interior. Moreover, as was pointed out in Q, there exist models with multiple necks. In 
fact the number of necks can be unbounded for a given equation of state. Some examples of double neck 
models with quite pronounced double necks were given in |q] . In the present paper we show how such models 
can be constructed. We also prove the statements made in [Q. 

The relation between necks and the optical geometry can be described in terms of the centrifugal potential 
iVi) which governs the null geodesic orbits in the stellar spacetime. A slight modification (V = Vi + Vm) of 
this potential is responsible for the propagation of axial gravitational wave modes. In this picture a bulge 
corresponds to a potential well in which trapped modes may exist. Models with multiple necks (and hence 
bulges) have potentials with multiple wells. The corresponding gravitational wave modes can therefore be 
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considerably more complicated than the single well modes. 

Our analysis relies on a recent dynamical systems approach by Nilsson and Uggla [^ Q who treated the 
full solution space of the SSS field equations, with focus on the case of linear and polytropic equation of 
state. Nilsson and Uggla give three different dynamical systems reformulations of the field equations, one 
for the linear case |q| and two complementary ones for the polytropic case m . Referring to the shape of the 
corresponding state spaces, we denote these formulations as the saddle systemR H, the loaf of bread system 
and the cube system [Q. 

Here, we will be using a variant of the loaf of bread formulation to prove the existence of optical geometries 
having an arbitrarily large number of multiple necks B. Surprisingly this phenomenon occurs for the 
extremely simple Zel'dovich stiff matter equation of state with an added constant, p = p — ps (we use the 
"s" subscript to denote evaluation at the stellar surface where p = 0). This equation of state is everywhere 
causal, the speed of sound being exactly equal to the speed of light. However, the necks are quite unobtrusive 
and apparently do not give rise to trapped gravitational wave modes. Therefore, for the sake of illustration, 
we use models with a first order phase transition in this paper. For simplicity we focus on models with 
two uniform density layers. Such models were recently considered by Lindblom ^ who studied their mass- 
radius curves. It turns out that the phase transition has a quite dramatic effect on the gravitational field 
in the stellar interior. In particular it is possible to achieve very pronounced necks if the phase transition is 
sufficiently strong g. 

We are leaving issues of stability aside in this paper. However, we may remark that the Zel'dovich 
equation of state is in fact unstable beyond the first neck as can be seen from its mass-radius curve. The 
stability of the two-zone uniform density models remains unclear. The mass-radius curves cannot be used 
since there is no well-defined adiabatic index given by the equation of state itself. Therefore a dynamical 
stability analysis must be employed to determine the stability. However, there is a technical difficulty here 
since the perturbation equation is not of the Sturm-Liouville type because of the discontinuity in the density 
function. Hence the standard theorem for the spectrum of eigenfrequencies does not apply in this case. 

2 The modified loaf of bread system with the inverse chemical 
potential as matter variable 

In the Schwarzschild radial gauge, the metric for a static spherically symmetric (SSS) spacetime can be 
written as 

ds2 ^ -e'^'^dt^ + e^^dr^ + r^{d9^ + sm^ 9 dcf^) (1) 

with v and A being functions of the radial variable only. In the case of an isentropic perfect fluid, described 
by a pressure p and an energy-density p constrained by an equation of state p ~ fip), the stress-energy 
tensor takes the form 

Tab = pUaUb+p{gab + UaUb), (2) 

where m° is the four-velocity of the fluid. The SSS symmetry implies that u°' = e~'^{d/dt)°' and that p 
and p are functions of the radius only. We shall restrict our attention to fluids satisfying the basic physical 
requirements that p, p and dp/dp are everywhere non-negative. We will also assume that the equation of 
state has the following asymptotic behaviours 

p/ p — > as p ^ (3) 

p/p -^ constant > as p — s- oo (4) 

The p — > behaviour imposes no restriction on stellar models, since it is necessary p| for the existence 
of solutions with flnitc radii. The p ^ oo behaviour could also be expressed as saying that the speed of 

^ The "saddle" here refers to the shape of the surface of vanishing pressure (which serves as a boundary for physical solutions) , 
rather than the shape of the full state space. 



sound Usound = \/ dp/dp should not go to zero (or oscillate) in the given limit. Furthermore, if p/ p were 
to vanish in the large pressure limit (so that p/p diverges), solutions with large central pressures would be 
Newtonian [p <C p) near the center, while, on the contrary, a large central pressure star is expected to be 
highly relativistic near its center. 

For future reference we note that condition (H) implies that the equation of state can be uniquely split 
in the useful form 

p = (ry - l)p + -0, (5) 

where 77 = 1 + limp_,oo p/p > 1 and ^/i is a function of p satisfying i/j/p —> a.s p —> 00. Furthermore eq. (0) 
implies that p/ip ^ as p — > 0. There are further physical conditions that can be imposed, specifically the 
dominant energy condition p < p, as well as the condition that the fluid should not become superluminal, 
which is saying that the velocity of sound Wsound should not exceed the value 1. However, these two conditions 
need not be stressed in this general discussion, but note that they independently imply r/ > 2. 

A global picture of the solution space of Einstein's equations Gab — f^-Tab for such a model can be obtained 
by formulating the equations as a dynamical system, preferably on a compact state space using dimcnsionless 
variables. In this work wc use the loaf of bread formulation of Nilsson and Uggla as our starting point. Later 
on we will modify that formulation to better suit our purposes. In the loaf of bread formulation Nilsson and 
Uggla use two geometric variables given by 

(6) 



1 + ri^^r 

g2A 



where v^r = dv /dr. In these variables, the quotient m/r is given by 

1-^^.-^^, (8) 

r K 

where m is the standard spherically symmetric mass function of Misner and Sharp |nQ|, which can be 
invariantly defined in terms of the Schwarzschild radius r by the relation 

1 ^VrVaT. (9) 

r 

For the matter degree of freedom, Nilsson and Uggla chose 

y^^- (10) 

p+p 

The independent variable x is defined up to translations by 

dx= ,/'. . (11) 



This gives the evolution equations 

S' = -yi^S + iP[l + 22/(1 - 2E)] (12) 

K' = 2y{l - K - 2P)K (13) 

y' = -yF{y)ll (14) 

where P := l — Yi^ — K, F{y) := y{p/p — dp/ dp) — y{ip/p — dip/dp) and a prime denotes differentiation with 
respect to x. This formulation is based on the assumption that the function F is strictly positive except at 



infinite pressure {y ^ rj ^) where it becomes zero. When the assumption holds, the boundaries of the state 
space are defined by the following invariant submanifolds 

y = (surface of zero pressure) (15) 

y — rj^ (surface of infinite pressure) (16) 

P = (i.e. a: = 1 - S2) (17) 

X = 0. (18) 

This clearly gives a compact state space, whose coordinate image in M"^ is shaped like a loaf of bread (see 
fig. H). There are no equilibrium points in the interior of the state space, implying that the topological orbit 
structure is completely determined by the equilibrium points on the boundary surfaces. The equilibrium 
points and their associated eigenvalues are given in [Q for the case of a polytropic equation of state. The 
evolution equation for the pressure being p — — pS, the right half of the state space defined by S > is the 
physical one where a positive pressure is non-decreasing. A physical solution is furthermore required to be 
regular at the center of the star, which implies that the center corresponds to a point on the "ridge" (S = 0, 
K = 1) oi the loaf. The surface of the star corresponds to a point on the arc y = 0, P = 0, S > 0, which 
serves as an attractor for all interior points of the state space. By the definition of the independent variable 
X, an orbit representing a physical star approaches a point on the ridge as x — > — oo and a point on the 
surface arc as a; ^ oo. By evaluating eq. (^ on the surface arc where K — 1 — S^, one finds M/R — mg/rs 
being given by 

which determines the matching to the exterior Schwarzschild geometry. An orbit for which Eg = has an 
infinite radius R, but could still have a finite mass A/, but here we do not address this issue in detail since 
we shall only be interested in models with finite radii. 

According to eq. (04) , the assumption of F being positive for finite pressures means that y is a monotone 
variable. However, although eqs. (H) and (0) ensure that F will be positive in the limit of small as well as 
large (but finite) pressures, a physically reasonable equation of state can have a soft intermediate regime 



where F drops below zero, a notable example being the Harrison- Wheeler equation of state |11 . For such 
equations of state, this formulation is only suited to describe stars for which the central pressure pc is less 
than the lowest pressure pi which gives F = 0. The reason is that every pressure pk for which F = defines 
an invariant submanifold y = yt- Hence for a star with p^ > Pi, the independent variable x will run off to 
infinity before the surface arc is reached. 

To remedy this situation we shall use a different matter variable, namely one with a simple relation to 
the manifestly monotone, bounded and dimensionless function 

w = e"'*, (20) 

where 

^-r^^- (21) 

Note that h has a simple relation to the gravitational potential, v = —h + constant. By its definition /i is a 
continuous function of the pressure whether or not the same holds true for the density p = f{p). Hence, by 
using a continuous function of h as matter variable, we will also be able to treat equations of state with a 
phase transition linipf p,. p < limp|pj p at one or several pressure values pt ■ 

The function w has a simple physical interpretation of being proportional to the inverse of the chemical 
potential, which often is defined as ||l2| 

,-'-^, (22) 



where rip is the particle number density. It then follows from the relation np = {dfi/dup)^^ ||l^ that the 
chemical potential may be reexpressed as /i = ^sW^^ ■ 

For the equations of state we consider, the function h has the asymptotic behaviours 

hp/p -^1 a.sp-^0 (23) 

h/ Inp -^ rj^^ as p -^ oo (24) 

Since by assumption p/ p — > as p — > 0, eq. (E3) shows that although the integrand of eq. ( po| ) may be 
divergent at p = 0, the integral itself is well-behaved. It follows that w has the compact range < w < 1. 
There are, however, two reasons why we shall not choose to use w itself as the matter variable for the 
dynamical system. Firstly, we want the matter variable to decrease with decreasing pressure and to become 
zero at p = 0. This will make our reformulation qualitatively more similar to the original formulation which 
uses y — p/{p -\- p)- The requirement can of course always be trivially met, e.g. by using 1 — u; rather than 
w. Secondly, if the matter variable is not chosen properly, the dynamical system will in general become 
irregular on the surface of infinite pressure. Since there will be equilibrium points on this surface, important 
for the solution structure, this would be an unwanted feature. 

Referring to our sought for matter variable as y, its evolution equation can be written as 

dY 
Y' = yw-—T.. (25) 

dw 

Combined with the evolution equations ( |l2| ) and (^3[) for the geometric variables, we see that the system 



will be linearizable if y and wdY/dw are C^ functions of Y. According to eq. (14) and (|25|), we have 



dy _ F _ ^ A ^26) 

dY w dY/dw w dY/dw ' 

where A ~ tp/p — dip /dp. Assuming that w dY/dw is nonzero at finite pressures, the above expression for 
dy/dY stays finite for finite pressures, but is clearly at risk of becoming divergent at infinite pressure, where 
y = ri~^ and A = w = 0. Hence we must choose Y such that the quotient 

^ = ^tlA (27) 

w dY/dw 

stays finite as p ^ oo. In order to avoid equilibrium points with zero eigenvalues, Q should preferably take a 
nonzero value in this limit. For simplicity we shall only consider the case when the leading behaviour of ip at 
high pressures is constant or polytropic, i.e. we shall assume that there exists a nonnegative finite number 
n such that 

c». (28) 
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(29) 

(30) 

in the same limit. It follows that we should not use a linear function of w as the matter variable unless 
n + 1 = rj. However, the choice 

y = l-^^ s = ^, (31) 

n + 1 

gives w dY/dw = — sw* — — s(l — Y) and hence a quotient Q which stays finite and nonzero as p — > oo. 
Finally settling for this choice, the evolution equations for the dynamical system take the form 

S' = -yKi: + iP[l + 2y{l - 2E)] (32) 

K' ^ 2y{l - K - 2P)K (33) 

y' = — f^y(i-ns, (34) 

n + 1 



where y should now be viewed as a function of Y . Strictly speaking this system is not a proper dynamical 
system if there is a phase transition. In that case y has a discontinuity which leads to discontinuities in the 
derivatives of all the state space variables. This is however not a serious problem as long as the numerics is 
done carefully near the phase transition. 

It may be noted that for an equation of state of the type p = {rj — l)p + fcp"/'"+^-', i.e. ip = fcp"/("+^\ 
our matter variable Y can be explicitly expressed as a function of p as 

In fact, this shows that Y is then proportional to the original matter variable y according to Y = rjy. This 
particular class of equations of state includes the incompressible case (ry = 1, n = 0), the polytropic case 
{f] — 1, n > 0), the relativistic polytrope {i] ^ n + 1, n > 0) as well as the linear case {n — 0, k — ps). More 
generally, for every equation of state that satisfies our general assumptions (Q) , (^ and ( p8| ) , the qualitative 
picture of the state space resembles that of the original formulation for the case of y being monotone. 
Quantitatively, the invariant boundary submanifolds y = and y — rj^^ are replaced by y = and Y = 1, 
respectively. We stress that there will be no additional submanifolds even if y is nonmonotonc. Moreover, 
the reformulation alters the eigenvalues and eigenvector components of the equilibrium points. The new 
eigenvalues of all possible equilibrium points are given in table H. 

As far as the numerical integration of the system given by eq. (p3), (p3) and (p3) is concerned, one 
encounters a difficulty when using an equation of state for which it is not possible to explicitly express y as 
a function of Y , since y has explicit occurrence in the left hand sides of the equations. In the case when the 
equation of state is given in parametrized form, 

p=p(A), p^p{\), (36) 

with p and p being monotone functions of the parameter A, one can easily get around the difficulty by viewing 
y as a function of A and utilizing A as an auxiliary dynamical variable with the evolution equation 

y. -,(!)-.. ,37, 

This gives two possibilities of calculating Y . Either one can integrate the three evolution equations for E, 
K and A and, given the result, calculate Y according to 

Y = l-e.J-^f'^dx], (38) 

y n+lJx, P + P J 

or else one can integrate the evolution equations for S, K, Y and A simultaneously. In the latter case eq. 
( pq ) plays the role of a propagated constraint, which offers the possibility of putting the numerics to a test 
by afterwards checking the accuracy to which the relation holds. 

For each equation of state, there is a one-parameter family of stars corresponding to the choice of central 
pressure Pc- This choice is the sole initial condition, specifying the starting point y = Y^ up on the ridge of 
the state space. However, since all points on the ridge are equilibrium points, the starting point of the actual 
integration has to be moved out a tiny distance from the ridge point in the eigendirection of the positive 
eigenvalue 2yc- This point is given by the following linearized solution near the center: 

E = 2(l + 2?;c)ee2«<^^ (39) 

K = l-12ycee^y'''= (40) 

Y^Y,-s{l-Y,){l + 2y,)ee^y^\ (41) 

where e is the trivial integration constant corresponding to the freedom to make translations in x. The 
corresponding linearization for the auxiliary variable A is 

l + 2t/c fdp\^ 2y,x f.^-. 

X^K-—^pJ-\ eey . (42) 



3 The surface of regular orbits 

The orbits with regular centers, i.e. the orbits starting out on the ridge E = 0, A' = 1 and leaving it 
according to the hnearized solution given by eqs. ( p^ ) - (pl|), defines a subset V of the state space S which 
we will here discuss in some detail. The fact that the interior of S contains no equilibrium points implies 
that all orbits begin and end at equilibrium points on the boundary dS. As a consequence, the interior of 
V is an embedded two-dimensional submanifold in the three-dimensional interior of 5'. To find out how this 
two-surface is embedded, it is clearly important to understand the way the boundary dV is situated in S. 
To this end we first need to understand the stability properties of the subset of equilibrium points in dS 
which also belongs to dV. This set of critical points is comprised of the line of stellar center points L2, the 
single point Pi on the surface of infinite pressure F = 1 as well as the line of stellar surface points Li (the 
part of the line of equilibrium points Li for which E < Emax for some E^ax < 4/5 which according to eq. 
( |l9| ) gives the minimal value of the tenuity R/M > 9/4). The eigenvalues of these points as well as the exact 
position of Pi can be found in table |l| and are shown in figure y. 

Let us first discuss the line of points L2, which is the part of dV from which the regular orbits start out. 
Setting Y = gives the point which connects L2 with Li and simply corresponds to the Minkowski solution. 
As can be read off from table m, this point is non- hyperbolic, i.e. all three eigenvalues vanish implying that 
the orbital structure near this point can not be obtained by a linearization of the system of ODEs. For 
< y < 1, points on L2 has one zero, one positive and one negative real eigenvalues. While the negative 
eigenvalue is associated with an eigenvector tangential to the boundary surface P = 0, the eigenvector of 
the positive eigenvalue, which generates the regular orbits and can be found in table g, is directed into the 
S > half of the interior of S and towards decreasing values of Y. However, when going to infinite central 
pressure by setting Y = 1, the Y component of the latter eigenvector vanishes and generates an orbit which 
necessarily is lying in the invariant submanifold Y — I and is ending at Pi. Since this orbit never reaches Li, 
it obviously corresponds to a stellar model which is infinite in extent. Indeed, since the form of the pressure 
function ip in the equation of state p — [rj — l)p + ^p is irrelevant as long as one stays on Y = 1, it gives a 
solution to the Einstein equations which can be identified with any regular center solution that is obtained 
for the scale invariant equation of state p = {rj — l)p. This follows from the fact that all regular center 
solutions for the scale invariant equation of state are equivalent as they simply are related by a rescaling 
of dimensionful variables, a change of units if one prefers, and hence they are all equivalent to the infinite 
central pressure limit solution. For this reason, we shall refer to the orbit in question as the RSI (regular 
scale invariant) orbit. Note that the RSI orbit is a part of dV, connecting L2 with Pi. 

As far as the equilibrium Pi is concerned, this point itself again corresponds to a solution for the scale 
invariant equation of state, namely the self-similar Tolman solution, which neither has a regular center nor a 
finite radius. The point Pi has two eigenvalues whose real parts are always negative and whose eigenvectors 
are tangential to the F = 1 surface. For rj < rj^, rj = rj,, and -q > -q,, with 77* = —2 + 8/\/7 ~ 1.024, these 
eigenvalues are real and distinct, real and equal respectively complex conjugate. In particular, this means 
that for rj > 77* the RSI orbit will spiral around Pi infinitely many times before reaching it. The third 
eigenvalue is always real and positive and the associated eigenvector, which is given in table 0, points into 
the interior of S. It generates an orbit corresponding to a solution which does not have a regular center but 
which does end on the stellar surface line Li. This orbit will be refered to as the skeleton orbit. Just like 
the RSI orbit, the skeleton orbit is part of dV , this time giving the connection between Pi and Li. It is 
important to understand that one may view the combined RSI and skeleton orbits, joined together at Pi, as 
a limiting orbit for the regular center orbits: by choosing a finite but sufficiently large central pressure, it is 
possible to obtain a regular orbit which follows this combined limiting orbit arbitrarily closely. In particular, 
let us assume that rj > rj,, and focus on the projection of finite central pressure orbits onto the Y = \ plane. 
Since the RSI part of the combined limiting orbit in this case is a spiral, there can be no upper limit on the 
number of times such an orbit may spiral around the point Pi, given that 77 > 77, . As we shall see, this is 
closely related to the possibility of constructing stellar models with an arbitrarily large number of necks in 
the optical geometry. 

We now finally turn to Li. Equilibrium points on Li with E > have one zero and two negative 
eigenvalues, Ai — — E and A2 = — PgE, which arc distinct miless the matter function F = y{tp /p — d'tjj/dp) 



takes the value 1 at zero pressure. Reexpressing the surface value of i^ as Fg = 1 — \iinp^o{p/p){dp/dp) 
reveals that < Fg < 1 and that Fg only depends on the leading behaviour of -0 in the limit of small 
pressures. We note that if the leading behaviour is polytropic, i.e. if there is some number n > such that 
^p-"/(n+i) _> const > as p ^ 0, F will evaluate to the number l/{n + 1) at zero pressure and hence will 
take the value 1 if n = 0. We note in passing that the polytropic type index n in general need not coincide 
with the index n defined by eq. (pq). Moreover, we remark that it is not true that F, = 1 iff the surface 
density ps > 0, since e.g. a leading behaviour of the perhaps unrealistic type ip ln(pi /p) — > const > as p -^ 
gives Fg = 1. If Fg < 1 so that Ai > A2, the eigenvector Vi belonging to Ai is directed in the pure S direction, 
while the eigenvector V2 of A2 is instead tangential to the state space boundary P = 0, its K component 
being —21] times its S component. From table it can be read off that the E component has the same 
sign as the Y component, which should be chosen negative for the eigenvector to point towards decreasing 
pressure. However, since all eigenvectors are here found to be tangential to state space boundaries which 
are invariant submanifolds, no interior orbit can enter a point on Li exactly along an eigenvector direction. 
Setting Fg = 1 so that Ai = A2, the two eigenvectors Vi and V2 are obviously still linearly independent which 
implies that every linear combination of these two vectors is also an eigenvector. In this case it is impossible 
to use the eigenvector structure to draw any conclusions about how the regular orbits enters Fi, except for 
the fact that they must enter tangentially to the plane spanned by Vi and V2. Some numerical studies using 
the equation of state p = (rj — l)p + fcp"/("+i) for which Fg — l/(n + 1) seem to indicate the following. 
For Fg equal or close to unity, the regular orbits in general do not tend to enter ii along neither Vi nor 
V2. However, with increasing Fg, the orbits tend to approach Fi tangentially to the P — surface, along 
the eigenvector direction provided by V2. This is to be expected, since the eigenvalue quotient A2/A1 = Fg 
decreases with Fg. Moreover, with increasing n, the orbits will tend to end at smaller values of E until 
eventually all orbits will end at the non-hyperbolic Minkowski point at E = where the local analysis breaks 
down. 

Now, an important point regarding the Fi part of Fi, is that its points need not be in one to one 
correspondence with the regular center orbits. Since the value of E on Fi by eq. (M) determines the 
compactness M/R, this is directly related to the fact that the mass-radius curve associated with a given 
equation of state does not in general define a monotone relation between M and R. In particular, the 
maximal surface value Emax (and hence the maximal compactness M/R) need not be the one obtained for 
the skeleton orbit, i.e. when going in the limit of infinite central pressure. Indeed, for many equations of state 
the mass-radius curve ends in a spiral, corresponding to the surface of regular orbits being rolled around 
its skeleton orbit boundary, always without crossing itself, before it is finally squashed into the Fi line. As 
we have seen, when 77 > 77*, the surface is already rolled around the skeleton orbit when leaving the Y = 1 
surface, so it seems likely that the high pressure regime of the mass-radius curve is closely related to the 
structure of the RSI orbit. 

4 Computing the optical geometry using the loaf of bread system 

In this section we will use the previous analysis to study the optical geometry]^ of static spherically symmetric 
models. The optical geometry may be defined as the spatial part of the conformally rescaled original 
geometry, with a specific conformal factor. The symmetries allow us to visualize essentially all important 
features as an embedded surface in a 3-dimensional Euclidean space. The geodesies of this surface then 
correspond to the null geodesies of the original spacetime which enable us to use our ordinary Euclidean 
intuition to investigate phenomena such as closed lightlike orbits and their stability. We now derive the 
equations that describe the embedding. 

We start from the Schwarzschild metric (||) and consider only the null geodesies, i.e. ds^ = using the 
fact that conformal rescalings of the metric leave them invariant. We may thus multiply both sides by the 
non-zero conformal factor e~^'^ . We also restrict attention, without loss of generality, to the equatorial plane 
9 = 7r/2 thereby obtaining the new metric 

= ds^ = -dt' + e^^^-^dr^ + 6-2^2 d02 (43) 
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Table 1: The equilibrium points with associated eigenvalues 



Eq. point 


Eigenvalue 


Eigenvector V ^ [V^ , V^ , V^] 
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[1,0,0] 
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[2(l + 2y),-12y,-s(l-y)(l + 2y)] 
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2 

(71+l)(77 + 2) 


[-2(?7 + 2){n(3 + 2?7 - 4), 16(77 + l)/3, -(r/ + 2)3 (/3s-i + 2.s + 77 + 2)] 



Table 2: Eigenvectors of some important eigenvalues. Definitions: (3 = rf + At] — 4, s = rj/{n + 1) 



It is clear that the spatial geodesic equations will be independent of time, and that the timelike part will 
have a trivial solution, so we drop the time part of the metric and turn our interest to only the spatial part 



df 



,2{X-u)^^2 _^ g-2.^2 



(44) 



This expression may be written in a standard form if we introduce two new variables given by dr* = e^^ "^'dr 
and r — e^'^r, yielding 



dl^ = d? 



■? + r\ 



(45) 



so that radial distances in the optical geometry are measured by dr, , while distances along a circle of constant 
r* are measured by fdcj). For this reason we refer to r* as the optical radius and f as the optical circumference 
variable. Even though r* is a monotonous function of the Schwarzschild radius r, the relationship need in no 
way be simple, so it is not possible to directly read off r from the optical geometry. For later reference we also 
point out that f has a simple connection to the centrifugal part of the effective potential, Vi = l{l + l)/f'^- 
We now wish to embed the optical geometry ( |45| ) in M^. To this end we write down the Euclidean metric in 
cylindrical coordinates. 



ds| = dz^ + dp^ + p^d(l)^ 



(46) 



The embedded surface will, at least locally, be given by some functional relation f{z,p) — 0, which, again 
locally, may be solved to give z = h(p), where h is some function. Thus dz = h^pdp, and the induced surface 



metric becomes 



dl' 



[hpY d~p' + p'd<j,'. (47) 



Comparing (B^) with (Es) we find the relations 

dh^ — drl — df^, p = r{r^), 

which defines a differential equation for h{p) given below. 



dh 



^^e-''^Je^^-{l-ly^rrf (48) 



dr 

In the interior of the star we may express this system in the state space variables of section 2. Thus extending 
the state space by the variables f and h we are able to solve the embedding equations simultaneous with the 
dynamical system. The new structure equations are 



r = yf(l-2I]), h = yry/K - (1 - 2E)2 (49) 

It is obvious that the center of the star (f — 0) still defines an invariant submanifold for this system, so we 
must again linearize the equations to find initial conditions at a finite distance from the center of the stellar 
object. The linearized solution is in this case 



^2 



eexp(2j/ca;) 



h = y^2 + yceexp{2ycx), (50) 

where we were forced to linearize f^ rather than f, since as x — > — oo, f^ is asymptotically a linear function 
of the original variables E, K and Y, whereas f itself is not. The solution thus obtained may be scaled to 
any given mass of the star by matching it to the appropriate exterior Schwarzschild solution at the surface. 
In the exterior we have e^" = 1 — 2^ and e^^ = (1 — 2C)~ , where C, = M/r, so the embedding equations 
( [fq ) take the form 



VT^2C' dr y (1-20^ 

From the fact that necks and bulges correspond to the zeros of f,r it is clear that C = 1/3 correspond to a 
neck in the exterior geometry iff i? < 3M. We may remark here that in the interior the necks and bulges 
correspond to the zeros of E — 1/2 as is easily seen from ([49|). 

5 Optical geometries with an unlimited number of necks 

Using the fact that an interior neck or bulge in the optical geometry corresponds to the state space orbit 
crossing the surface E = 1/2, we now prove what was claimed in B, namely the existence of causal models 
for which the number of interior necks can be arbitrarily large. To see this, we begin by reading off from 
table n^ that the position of the Tolman equilibrium point Pi on the Y = 1 plane depends on the equation 
of state parameter rj. In particular, since the E coordinate of the point is 2/{r] + 2), Pi can be put on the 
E = 1/2 surface by setting rj = 2, which puts it at E = K = 1/2 and gives a speed of sound that approaches 
that of light in the limit of infinite pressure. The assertion now follows almost immediately from our previous 
discussion of the surface of regular orbits. Since ry = 2 exceeds the critical value 77* = —2 + 8/a/7 « 1.02 
above which the RSI orbit spirals into Pi, an arbitrarily large number of E = 1/2 crossings for a regular 
and finite central density orbit can be obtained by choosing a sufficiently large (but finite) central pressure, 
making the orbit follow the combined RSI/skeleton limiting orbit sufficiently closely. Although this proof 
does not depend on the particular form of the pressure function ^ in the equation of state, an arbitrarily large 
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Figure 1: The state space of the loaf of bread system for the stiff equation of state p = p+Ps- The equihbrium 
points, indicated as bold dots in the figure, are denoted as in table |l|. The figure also displays the combined 
limiting orbit described in section 0, corresponding to the limit of infinite central pressure within the family 
of regular center orbits. The RSI part of this combined orbit, which completely lies in the boundary y = 1, 
starts out from the Y = 1 point on L2 and spirals into Pi. However, from the figure it is not at all clear 
that the RSI orbit actually spirals into Pi, which is explained by the fact that the absolute value of the 
quotient q between the imaginary and the negative real part of the relevant complex conjugate eigenvalues 
is not very large, the precise value being q ~ ^/3. This implies that the distance to Pi decreases by a factor 
g-27r/g r^ 0.026 for every turn around Pi. The skeleton part of the combined orbit, around which the surface 
of regular orbits is rolled, is for this particular equation of state confined to the plane S = 1/2. 

number of necks may not always be obtained for models with finite radii. However, since linear equations of 
state p — (rj — l)p + ps (ip = ps = const > 0) are known to always give regular solutions with finite radii, we 
may simply choose the stiff equation of state p = p + ps to avoid the possibility of ending up with models 
that are infinite in extent at large central pressures. We also remark that for the stiff equation of state, the 
eigenvector generating the skeleton solution has a vanishing E component, as can be read off from table 0. 
In fact, the skeleton solution is in this case an exact solution which is entirely confined to the S — 1/2 plane, 
namely the solution E = 1/2, K = |(?/ + 1)~^ [0|. This is likely to imply that the number of necks grows 
faster with the central pressure than it does for more general forms of ip, since a nonvanishing E component 
of course drives the skeleton orbit away from E = 1/2. The combined RSI/skeleton orbit for the Zel'dovich 
equation of state is shown in fig. H. 

6 Some multiple neck optical geometries 

Here we investigate the interior solution for different equations of state. Our present purpose is to show 
that the interior optical geometry will look quite different for the different cases. We showed in an earlier 
paper |^] that the geometry may develop additional "necks" and "bulges" corresponding to unstable and 
stable closed lightlike orbits respectively. As shown by Chandrasekhar and Ferrari ||lj] the axial modes of 
nonradial metric perturbations of SSS perfect fluids are (to first order) not coupled to a perturbation of the 
matter. Such metric perturbations may therefore be interpreted as pure gravitational waves. By separation 
of variables the equations governing these modes can be reduced to the 1-dimensional Schrodinger equation 



dr: 



-Zi + VZi = a'Zi 



(52) 
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where the total potential V = Vi + Vm is dominated by Vi = l{l + l)r^^, I = 2,3, . . . which is simply the 
effective potential for photons with angular momentum L^ ~ 1(1 + 1). This, together with the fact that the 
correction Vm to the total potential is small compared to Vi in most cases, shows that the optical geometry 
is closely related to the trapped modes of gravitational radiation since the maxima and minima of Vi clearly 
correspond to the necks and bulges of the optical geometry respectively. 

We turn, now, to the analysis of the optical geometry of some different models. In the case of a double 
layer uniform density model the equation of state is parametrized by the constant densities p-^- and p_ 
of the interior and exterior zones respectively, as well as the transition pressure pt where the density is 
discontinuous. In the single layer uniform density model considered by Abramowicz et. al. H the only 
parameter is the density po, and the linear model is parametrized by the central and surface density pc and 
ps respectively. Given these parameters, each stellar model is obtained by specifying the central pressure Pc- 

Below we display and compare the optical geometry for these models. The parameter values for the 
two-zone models are chosen to obtain three distinct examples of optical geometries with pronounced double 
necks, corresponding to two deep potential wells as shown in [^. They were chosen to have the same value 
of the tenuity a :— R/M = 2.2749, as well as the same mass. This condition is visualized by letting the 
exterior geometry be identical in each triplet of figures. The chosen parameters are shown in table ^, where 
we have chosen to use cgs- units and assigned a mass of 1.5Mq to each model. We see that the central density 
required to produce the double neck structure is very large. We could have chosen to have a more realistic 
central density, but then the total mass of the corresponding stellar object would have been of the order of 
several hundred solar masses. 

In figure we compare a two-zone model with a uniform density model, and a model with a linear equation 
of state. The parameters of the two-zone model were chosen so that the bulges should have approximately 
the same size corresponding to two equally shaped potential wells. The density of the uniform model were 
then chosen to give the same tenuity as this model. The parameters of the linear model were chosen so that 
the single neck should be reasonably visible. The tenuity of this model is a = 2.81 which is close to the limit 
cimin = 2.749 for this type of equations of state. The lower part of the double layer model is of course, as 
in the uniform case, a sphere. We thus expect that the modes of the inner potential well, granted that the 
inner neck is thin, should resemble those of the uniform model. 

In figure we look at three different two-zone models. The first is the one already shown in figure The 
two other models are chosen to display the possibility to have either a large upper bulge or a large lower one. 

It would be interesting to know if the two-zone models are stable against small radial perturbations, but 
this question turns out to be nontrivial. Of course, a star with constant density is always stable, however 
since this configuration is very unphysical (the speed of sound diverges) one usually views the equation of 
state p = constant as an approximation to the more physical p = pa + pi{p), where po is a constant and 
pi{p) is a function of the pressure satisfying pi{p) <C po for all p < pc. One may then impose reasonable 
conditions on this equation of state, for instance causality (usound < 1), and the dominant energy condition 
{p "^ p)- It seems however that you cannot find any reasonable functions pi{p) such that these conditions 
hold and the star is ultrarelativistic. We may note that none of our two-zone models respects the dominant 
energy condition. Ignoring these physical difficulties and viewing the models as toys, we may proceed with 
the stability analysis. Firstly we may remark that you cannot use the usual analysis of the mass to radius 
curve to establish stability properties for any model with constant density since the adiabatic index is not 
defined |l^ (it is determined by the small function pi(j)) which is arbitrary here). We are therefore forced 
to try a dynamical approach, but here the analysis is complicated by the discontinuity in the perturbation 
equation so that the differential equation is no longer of the standard Sturm-Liouville type. Therefore, the 
standard theorems of Sturm-Liouville theory cannot be directly applied. However, since the discontinuous 
density function can be arbitrarily well approximated by a continuous function one might expect that nothing 
too extraordinary happens. Preliminary results using a variational principle indicate that this is indeed the 
case. The stability of more realistic ultracompact (in the sense R < 3M) stellar configurations with a 
phase transition has previously been studied by Iyer, Vishvehvara and Dhurandhar p6[ . They conclude that 
there exist stable causal ultracompact models, however not with any of the high density equations of state 
proposed in the litterature at that time. 
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Two-zone models 


p+{g/cm^) 


P~{alcvr?) 


Ptig/cms^) 


Pcig/cms^) 


n/R 


Peanut 


2.74 X W^ 


5.48 X IQis 


3.52 X 10^8 


5.67 X lO-*! 


0.014 


Bigup 


5.81 X 10^^ 


5.51 X 10^5 


2.90 X 10^8 


5.22 X 10''! 


0.010 


Bigdown 


2.80 X W^ 


5.48 X 10^5 


3.63 X 1038 


2.90 X 10*2 


0.014 



Other models 


Pc{g/cm?) 


Ps(ff/cm3) 


Pc{g/cms^) 


Uniform density 


5.56 X 10^5 




7.59 X 10^7 


Linear 


9.00 X 10^6 


5.56 X lOi'^ 


7.59 X 10^7 



Table 3: Parameter values for the models considered. Here we have chosen the mass of the object to be 
1.5M0 which means that the radius will be about 6.2 km for the linear model, and 5.0 km for the others. 
We see that the term ultracompact is indeed relevant since these objects are denser than normal neutron 
stars are believed to be. We see that the two-zone models consist of very dense, small cores (rt refers to the 
radius where the phase transition takes place), and envelopes of more realistic densities. 

7 Concluding remarks 

The two-zone models discussed in this paper could be considered as toy models for stars having a first order 
phase transition in their interiors. Our results indicate that a phase transition is an important mechanism 
for producing a pronounced neck in the optical geometry of a stellar model. In this picture, the outer neck 
at r = 3M (if it exists) is associated with the vacuum/matter phase transition at the surface of the star. 
However, this is not the only mechanism for producing necks. Another example is provided by the stiff linear 
equation of state. In that case there is an increasingly long sequence of shallow necks which is associated 
with the high density spiral equilibrium point Pi B . 

The existence of multiple neck optical geometries is associated with gravitational perturbation potentials 
with multiple wells [pi . It would therefore be of interest to calculate the corresponding quasi-normal modes 
along the lines of calculations done for single well potentials in ||l| (uniform density models) and in p^ 
(polytropic models). 
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Figure 3: The optical geometry of the "peanut" model compared with the "bigup" and "bigdown" models. 
All models have the same mass and hence identical exterior geometry. 
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